Extraorinary Hall effect in hybrid ferromagnetic/superconductor (F/S) bilayer 
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Extraordinary Hall effect (EHE) in bilayer F/S(N) was investigated theoretically. The conduc- 
tivity tensor <t Q( 3 is calculated in the Kubo formalism with Green functions found as the solutions 
of the Gorkov equations. We considered diffuse transport in the ferromagnetic layer, taking into 
account as a main mechanism of electron resistivity s-d scattering. In this model Gorkov equations 
for s-electrons in the ferromagnetic layer remain linear and are solved easily. It is shown that Hall 
field E H for both F/S and F/N contacts are step-functions of the coordinate perpendicular to the 
planes of the layers and have zero value in S(N) layer. The Andreev reflection increases the value 
of Hall constant R s for F/S case. The value of the Hall constant is R F J S = i^ ulk (cr T + a 1 ) 2 /4<r T a 1 , 
where cr T and r are conductivities of electrons with up and down spins, and 7?^ ulk is the Hall con- 
stant in the bulk ferromagnetic metal. In fact, R^J S coincides with EHE constant of the bilayer of 
two ferromagnetic metals with equal thickness and opposite directions of their magnetizations. So 
we can make a conclusion, that the ideal interface between ferromagnetic metal and superconductor 
may be considered like a mirror with inversion in spin space. 
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The perpendicular spin-dependent transport in hybrid 
ferromagnetic/superconductor (F/S) and F1/F2/S struc- 
tures has been previously investigated in Ref. Jl], ||, [| . It 
was shown that the Andreev reflection j| at F/S inter- 
face causes a mixing of up and down spin channels and 
simultaneously a large spin accumulation on F/S inter- 
face arises. The total influence of both effects increases 
the resistance of the system at T = comparing to its 
value when the superconductor is in the normal state. In 
the case ofF^/F^S these effects strongly influence the gi- 
ant magnetoresistance (GMR) of the F1/F2 bilayer. For 
certain values of the parameters, the GMR could be even 
completely suppressed |§, || . 

Another spin-dependent transport effect, which has 
not been considered so far in F/S sandwiches, is the ex- 
traordinary Hall effect (EHE) §, §. The current j is 
assumed to flow perpendicular to the interface, the mag- 
netization M in the F-layer is in plane, and the Hall 
electrical field E H arises in plane in the direction per- 
pendicular to M . In such geometry, the Andreev reflec- 
tion is expected to play a more complicated role than on 
the current-perpendicular-to-plane (CPP) GMR because 
EHE combines current-in-plane (CIP) and CPP features. 

It is the purpose of this letter to investigate theoret- 
ically the EHE in F/S(N) bilayers in the diffuse regime 
in both situations where the non-magnetic layer is in its 
superconducting or normal states. In the present model, 
we assume that conductivities for spin-up and spin-down 
channels are different and the contribution of F/S inter- 
face resistance to the total resistance is small. This is 



the situation for which the CPP-GMR for F1/F2/S is 
destroyed by the Andreev reflection and spin accumula- 
tion 



I. GENERAL DEFINITIONS 

For geometry used in the model, the following relation 
proposed by Smith and Sears can be written 



Eg = j z (R H y + R s M y ), 



(1) 
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Ro is the normal Hall effect coefficient whereas R s rep- 
resents the EHE coefficient. H y and M y are respectively 
the external applied field and magnetization. The bias 
voltage is applied between the planes with coordinates 
z = —a and z — b, and z = is the position of interface. 

As it was shown in [|, || extraordinary Hall effect has 
two different origins — skew-scattering and side-jump. 
As a first approach to EHE in these structures, we will 
investigate only the first one, taking into account elastic 
scattering on impurities. The approach that we use for 
calculation of the Hall coefficient for layered system has 
been described in Following the same method, we 
calculate the components of the conductivity tensor using 
Kubo formula 

he 2 

<r al 3 = ^yTr [v a (G+ - G~) v (G+ - G~)] (2) 

which includes simple "bubble" diagonal conductiv- 
ity and vertex corrections for off-diagonal components, 
which are responsible for the extraordinary Hall effect. 
In (J2j> v a , V/3 are velocity components, indices a, (3 take 
the values x, y, z, G + ( - ' are retarded and advanced 
Green functions, is the volume of the system. 
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For the system under consideration, which is homoge- 
neous in xy-plane and inhomogeneous in z direction, it is 
convenient to use (re, ^-representation, where hn is the 
in-plane electron momentum. The Green function G is 
defined by equations: 



G = G° + G°H so G° 
G° = G e ff + GcffTGcs 



(3) 



where G° is the Green function of the system in the ab- 
sence of spin-orbit interaction, G e g is the effective Green 
function, diagonal on in-plane vector re, calculated in the 
coherent potential approximation (CPA) pi. T is a scat- 
tering matrix and H so is the spin-orbit interaction. If to 
adopt for the ferromagnetic layer the model of totally 
disordered binary alloy A c Bi_ c (c is a concentration of 
the alloy's component A) we can write the expressions 
for T-matrix in singe-site approximation and for matrix 
elements of H so in explicit form: 



S a u(n,z) - S c 



1 - [5°v(n,z) - g° + {z,z) N 

n 

v(n,z)i\M y [kxk'] y = ±J2HZ>( 



(4) 

-V t a+ ,(n,z) 
(5) 



In (U) and (||) 8" — e A — e B is scattering parameter, e A 
and e B are the band centers for A and B components 
depending on spin a, v{n 1 z) is projection operator: 



v(n, z) = as(n, z)c — ayt(n, z)(l — c) 



(6) 



1, if the site (p n , z) is occupied by 
a a (n, z) — ^ a atom (7) 
0, in the opposite case 

A = Xa — Ab is parameter of spin-orbit scattering which 
is non-zero only in the ferromagnetic layer. E CT+ is the 
coherent potential which is the solution of self-consistent 
equation (T CT+ (z)) = 0, where (. . . ) means averaging 
on impurities distribution, Q a + {z, z) = G^. + (z, z). 

In (J|) the z-component of electron momentum vector k 
in [k x k'] y is the antisymmetric gradient operator: k z — 

i (Vz - V 2 ). 

It is important to note that from definitions (||) and (j^) 
the first non-zero contribution into the vertex correc- 
tion of formula (|^) linear on H so is that containing 
(tl+(n,z)H™,(n>,z>))~S nn ,6(z-z>). 

In the adopted geometry the system of equations for 



Hall fields can be written as follows: 

j x (z) = J a xx {z,z')E*(z')dz' + 

+ Jj a xz (z,z',z")E z (z')dz'dz" = (8a) 
jz(z)= J a zz (z,z')E z (z')dz / + 

+ [[ a zx {z,z',z")E*(z')dz'dz" (8b) 



We consider H s ° /S like a small parameter of the theory. 
So a zx <C a zz and the second term in equation (8b) can 
be omitted. The off-diagonal component of conductivity 
has a three-point character. The additional coordinate 
z" represents the scattering plane. 



II. MODEL 

We consider a bilayer of the type F/S, where F is fer- 
romagnetic layer, S is a superconducting layer. A sim- 
ple two band (spin up and down) free electron model 
is adopted for this calculation. The Hamiltonian of the 
system is therefore written as: 

H = H F + H s (9a) 



H F = £ 

<r=+(4.) ,-(T) reF 



2m 



e F + sign(cr)e e 



x CWCW + 7 s dW «*(r)^(r) + ft.c) 



d 3 r (9b) 



H s = 



res 



T - 

^ \ 2m 



+ (A(r)^f* (r)Vf {r) + h.c 



d 3 r (9c) 



where £ c 



T2 12 

P F ~Pf 
2m 



is the exchange energy, ep, p F 



are respectively the Fermi energy and momentum. The 
second term in (|I]) describes the scattering of quasi free 
s-electrons into almost localized d-states. In bulk ferro- 
magnetic metals d-states may give contribute to the cur- 
rent fill . However in the present situation, we consider 
that there are no d-states in the superconductor. There- 
fore d-electrons are completely reflected on F/S interface 
and do not contribute to the current. On other respects, 
s-d scattering in ferromagnetic dirty d-mctal alloys re- 
mains the most important mechanism of s-electrons scat- 
tering @. In the case under consideration, we take it 
into account and consider that the random s-d scattering 
potential "f s d is much smaller then ep. We further calcu- 
late the mean free path in the Born approximation. A 
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is the order parameter in the superconductor. Now the 
system of Gorkov equations |12| for the normal G e g and 
anomalous F c g Green functions can be written: 



ll<iG\\{z,z) 



h 2 ( d 2 \ 

2^[d Z - 2 - H ) +£F + £ ^-^ dG ^ 

x GlUz, z') + AF%(z, z') = 5(z - z>) (10a) 



Perpendicular transport in bilayer F/S was investi- 
gated in H so we remind here only the result: E[l\ = 
E^l2, E\ and E[ are effective electrical fields acting on 
the carriers with spin up, down, e] = ,, , 2V ^ a } , , 

o(ll+i2) + f (1(2 ' 

j z ~ iqZi. Off-diagonal components of conductivity 
in (||) as well as diagonal one has two contributions - 
normal and anomalous. For example, one has the form: 



A*GU(z,z') 



d 2 



■ Ep — e e 



2m \dz 2 

+7 2 sd G^(z,z) F£(z,z') = (10b) 



We ommited index "eff" for the brevity. The terms £ e 
and 7 2 G 



2-oTT(ll) 



dil 



are different from zero and A = if z 
belongs to F-layer and vice versa in S-layers. The sys- 
tem ( |l0| ) is written for spin f . For spin J., £ ox has to be 

changed to -£ ox and G^ T) to G^ (U) . The main dif- 



ference between system ( |10D and usually employed equa- 
tion for F/S structures (see for example [13]) is that we 
took into account s-d scattering. Moreover, we consid- 
ered this s-d scattering as the main mechanism determin- 
ing the mean free path of s-electrons p"T[ . The function 
Gj%(z,z) may be considered like a constant in z-space 
if it is averaged over the short wave l eng th (h/pp) oscil- 
lations and the system of equation (nfl) may be solved 
analytically. Further, we set e cx = for s-electron. The 
explicit expression for Green functions are: 



Gll + (-oo < z,z' < 0) = 



p ifci \z— z'\ 

2iki 



F$ + (-oo < z,z' < 0) = 



fci 



G 



T1 - 



(— oo < z < < z < oo) = 



F$ + (-oo < z < < z < oo) 



G^ + (0 < z,z' < oo) 



i(ki + fc s ) 



(11a) 
(lib) 
(11c) 
(lid) 



4ifc, 



gik s z^—ik*z | ^ikslz — z 



^ik s z'^—ik*z ik*\z — z 



(lie) 



F$ + (0 < z,z' < oo) 



4ifc, 



gik s z^—ik*z | gik 3 \z— z | 



e ik s z' e -ik* 



ik* \z— z | 



(llf) 



where fci 2 = \ k 2 F - x 2 + ij^ = c 1:2 + idi*'> fcs = 



\Jk 2 F — n 2 + i\A 2 \ = c s + id s , li 2 are mean free paths 
for up and down spins. 



{a™™(z,z',z")U~L(v x „Gt(z,z") 
(T+(z") + H*°(z"))^, G+(z",z>)v z «,G-,(z',z") 

(T-(z")+^"(/'))^G-(z",z)) (12) 



where L means linear on H so part of this expression. 
Since T+„, = T+„ and H™, = -H%„, term © is 
proportional to (2i7? so (z)ImT + (z)) as well as all other 
contributions to a xz . In the Born approximation, we can 
write down for disordered binary system: 



(2iH so {z)lmT + (z)^ = -XM y c(l - c)(l - 2c) {S a f x 
x Im£r + (z,z)[fc x k'] y (13) 



from which it is easy to see that this term is proportional 
to (1 - 2c) /la- 

The unknown Hall field has to be found as a solution 
of the integral equation (|8a|). To solve it we take as a 
probe function for E H (z) the step function, taking value 



E 



HUD 



1 inside the ferromagnetic layer and in the 

superconductor for each direction of the electron's spin 
], {. In this case the equation (8b) may be rewritten as 
a system of two equations: 



^Iwd 1 — > — V'* 



2d x z 



-2d 2 a\ 



fclc 



dx - 



c 2 d 



(--"' • (1 - C 

3 
F 



-2d,b\ 



- J dx; ~ a < z < (14a) 



-2d\{z J ra) J2d\z^—2d2a^ 
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?Hl, 



-a«*.- (i _ e -a*«) dx . 



-2d, 



1-e 



-2d 2 a\ 



e -2d s (z+b) e . 



F 

2d B (z-b) 



dx = 



-2d\a Jld^a 



dx; 0< z <b (14b) 



where i?!? ulk is the Hall coefficient and crT(l) are conduc- 
tivities of up and down spin channels for the bulk ferro- 
magnet. For spin down we have to change T^-kLj 1^2. 

Solving the system of equations ( |l4a|) , ( |l4b|) we 
found that indeed the solution for the Hall field E H (z) 
in th e for m of t he step-like function satisfies the sys- 
tem (|14a|) , (14b) and consequently the integral equa- 
tion Q8a|) for any z. The results of the Hall fields are: 



£ft =E?i= 0, E» = E»%, E[ 



? H t _l E H 1 



fl bulkgT M ^ cr + o-- _ go the extraor dinary Hall coeffi- 
cient for F/S bilayer is: 



III. F/N BILAYER 

Now we recalculate the Hall coefficient when the non- 
magnetic layer is in the normal state. In this case, all 
Green functions are diagonal in the spin space: 

pikl\z-z'\ 

G[l + (-°° < 2, 2' < 0) = — — (16a) 



G 1 2 + (-oo < z < < z' < 00) = 



2iki 



(16b) 



i(ki + fc 3 ) 

tt , e ik 3 \z-z'\ 

Gil + (0 < z, z' < 00) = (16c) 

2lK3 



The system of equations for the Hall fields assumed to 
be step-functions can be written in the form: 



H ^ IX 3 / e -2di (z+a) g2di2 



1 1 1 klc\ 2 2 



l/2£f T Z 3 | 



fc F C v ' 



4 ./ fe|c 



R b s mk E{hM y 

,-2c2i(z+a) 2diz 



X 1 







1/2E?% [ -^e~ 2 <Wl 
J ktc 



^ dx; -a < z < (18a) 



e- 2dlQ ) d) 



E^h / 73 — I 1 1 dx = 



k pC \ 2 2 



-2d 3 z 



4 ./ fc|,C 



(l - e ~ 2<iia ) dx; 0<z<b (18b) 



of this system gives us E 2 = E 2 = 



For spin down we have to change T^->li 1 ^ 2. Solution 



H] _ VU*l + *L)M y ^ ulk (19a) 



4(0/3 + ^1 + ^1/3)" 



H1 ^ V^+^M, buU ., 
4(oZ 3 + bl 2 + %hh) ' 



(19b) 



where V is the total voltage drop across the F/N bilayer, 
and Hall coefficient is: 



(h + l 2 ) 2al 3 + b(h + h) + §Z 3 (Zi + h) obu 

it. 



2 a / 3 (/i + Z 2 ) + 2&W2 + |W2^3 



nlk 



(20) 



With the same method it is easy to show that the 
Hall constant for a spin-valve bilayer in antiparallel 
magnetic configuration F^/F^ is equal to Rf 1 ^ 2 — 

(^ T +^) 2 Dbulk 



40-To-i 



-R D . 



2kp 

where fc 3 = \jk 2 F — x 2 + i— — =03+ id$, 1$ is the mean 

V h 
free path in the normal metal. 

The following expressions are then obtained: 



4 _ fT al 3 + bh + ±hl 3 



E[ = E[ 



a/3 + bh + %hh 



E 



T _ 



Vh 

2{ah + bh + Uih) 



(17a) 



(17b) 



IV. CONCLUSIONS 



As can be see from (|l^) and (pOj), the relative change 
of the Hall coefficient in the presence of superconducting 
contact is equal: 

RP S -R^ _{o ] -o l ) 2 



and 



^buik 4crT(Ti 
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Simultaneously it is shown that resistivity p F / s of F/S 
bilayer is equal p F ^ s = ^ l — ^ 2 , where p\ and p2 are re- 
sistivities of up and down spin channels, this resistivity 
as well as R s coincides correspondingly with resistivity 
and the Hall constant of the bilayer of two ferromagnetic 
metals with opposite direction of magnetizations. There- 
fore, we conclude that an ideal interface between a ferro- 
magnetic metal and a superconductor can be considered 
like a quantum mirror with inversion in spin-space. Of 
course, the roughness of the interface may spoil the mir- 
ror image. In addition we considered the case where the 
spin-diffusion length was much larger than the thickness 
of the ferromagnetic layer. The influence of spin- flip pro- 
cesses on Hall effect will be considered in a forth coming 
paper. 

Experimental investigation of EHE in the situation, 



close to one described in the letter, may be done if to 
use as a ferromagnetic layer the alloy CuNi with rela- 
tively small exchange splitting and high resistance. The 
thickness i CuNl of CuNi layer has to be in the interval 
l e i •C t CuNl <c Z s d, where Z s d is the spin-diffusion length 
and l e \ is the longest mean free path of the electron in 
bulk CuNi. 
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